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We investigate low-energy fluctuations in the real kink crystal phase of dense quark matter within 
the Nambu-Jona-Lasinio model. The modulated chiral condensate breaks both the translational 
symmetry and chiral symmetry spontaneously, which leads to the appearance of phonons and pi¬ 
ons that are dominant degrees of freedom in the infrared. Using the Ginzburg-Landau expansion 
near the Lifshitz point, we derive elastic free energies for phonons and pions in dependence on the 
temperature and chemical potential. We show that the one-dimensional modulation is destroyed 
by thermal fluctuations of phonons at nonzero temperature and compute the exponent that char¬ 
acterizes the anisotropic algebraic decay of quasicondensate correlations at long distance. We also 
estimate finite-volume effects on the stability of the real kink crystal and briefly discuss the possi¬ 
bility of its existence in neutron stars. 


I. INTRODUCTION 

Revealing the nature of QCD at finite temperature and 
baryon density has been one of the most important chal¬ 
lenges in high-energy physics. At low temperature and 
low baryon chemical potential, the ground state of QCD 
is characterized by chiral symmetry breaking and color 
confinement. Under extreme conditions, QCD is known 
to exhibit novel phenomena, such as deconfinement of 
color degrees of freedom at high temperature and color 
superconductivity at high baryon density. These areas 
are extensively investigated; see [U3] for reviews. 

Inhomogeneous phases in QCD have been studied for a 
long time. Pivotal examples include p-wave pion conden¬ 
sation in nuclear matter [4] , a chiral density wave at large 
Nc [5“9] and crystalline color superconductivity at high 
density [10, 11], in close analogy to the Fulde-FerrelU 
Larkin-Ovchinnikov (FFLO) phases of superconductiv¬ 
ity [12, 13[ . Recently, an exotic phase of dense quark mat¬ 
ter in which the chiral condensate forms a spatial soliton 
lattice has attracted considerable attention [14-19], as re¬ 
viewed in [20]. Initially, a theoretical breakthrough was 
made in the study of inhomogeneous chiral condensation 
in (l+l)-dimensional field theories [21]. A general com¬ 
plex kink solution to the Ginzburg-Landau (GL) equa¬ 
tion for the chiral Gross-Neveu model was then discov¬ 
ered [22, 23] and its relevance in the phase diagram of the 
(chiral) Gross-Neveu model was elucidated [24] . The is¬ 
sue of inhomogeneous condensation in (3+l)-dimensional 
theories was revisited by Nickel [15] who introduced a 
novel technique for applying the general solution of the 
Gross-Neveu model to the Nambu-Jona-Lasinio (NJL) 
model and thereby pointed out the emergence of a one- 
dimensionally periodic chiral condensate in the vicinity 
of the first-order chiral transition line that appears when 
limiting to a homogeneous condensate. Gonsidering that 
higher-dimensional modulations are energetically disfa¬ 
vored against one-dimensional modulations [17, 18], there 
is a high chance that one-dimensionally modulated chiral 
condensates may indeed be the ground state of QGD in 
a certain range of temperature and chemical potential. 


So far almost all studies on the inhomogeneous chi¬ 
ral condensation have been limited to the mean-field ap¬ 
proximation and inclusion of collective fluctuations of the 
order parameter is an urgent issue. It has been well 
known since Landau’s seminal work in the 1930s that in 
three dimensions a one-dimensional spatial order is un¬ 
stable against thermal fluctuations. This is known as the 
Landau-Peierls theorem [25, 26] . The destabilized long- 
range order leaves its imprint in the power-law decay of 
the order parameter, characteristic of a phase with quasi- 
long-range order. This is reminiscent of two-dimensional 
systems with continuous symmetry, in which a long- 
range order is prohibited [27, 28] and the Berezinskii- 
Kosterlitz-Thouless phase emerges at low temperature 
[29, 30] . 

In this paper, we elucidate various properties of gap¬ 
less excitations on a spatially modulated chiral conden¬ 
sate in the NJL model. Our analysis is based on the 
6th-order GL expansion of the NJL model near the crit¬ 
ical point, for which an analytical solution to the GL 
equation, called the real kink crystal [22-24], is available. 
Because the kink crystal breaks both the translational 
symmetry along one axis and chiral symmetry, we en¬ 
counter smectic phonons in addition to ordinary pions. 
We will derive the elastic free energy of these modes and 
reveal how they modify the mean-field picture of this 
phase. While the analysis presented here closely paral¬ 
lels preceding works about fluctuations in liquid crystals 
[31-33] and in the FFLO phases of fermionic superfluids 
[34-39], our analysis has a specific focus on quantitative 
understanding of those gapless modes in the context of 
dense QGD. 

The outline of this paper is as follows. In Sec. II, we 
develop a general argument for low-energy fluctuations 
over a one-dimensionally modulated order parameter. In 
Sec. Ill, the GL-expanded NJL model is analyzed. We 
first derive the free energy for the translational phonon 
mode and argue that thermal fluctuations of phonons 
wash out the modulated condensate and lead to a phase 
with quasi-long-range order. Next we derive the free en¬ 
ergy of pions. Finally, Sec. IV is devoted to concluding 
remarks. Some technicalities are summarized in the Ap- 
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pendixes. 


II. SYMMETRY CONSIDERATION 

When a chiral condensate is modnlated along one di¬ 
mension, the phonon mode (u) appears as the Nambn- 
Goldstone (NG) mode of translational symmetry break¬ 
ing, in addition to pions (tt). In the real kink crys¬ 
tal phase discnssed later, the vectorial isospin symme¬ 
try SU(2)y is nnbroken and there is no mixing between 
phonons and pions. Therefore, the symmetry breaking 
pattern reads^ 

XI SO(3) — >■ [R^ XI SO(2)] X [discrete symmetry] 

( 1 ) 

in addition to the chiral symmetry breaking pattern, 
SU(2)i^ X SU(2)l ^ SU(2)v. Here R^* and SO(d) denote 
the d-dimensional translational and rotational symme¬ 
try groups, respectively. The discrete symmetry includes 
a discrete translational symmetry along the modulated 
direction of the condensate as a subgroup, which is a 
remnant of the translation R in this direction. Other el¬ 
ements of the discrete group depend on the shape of the 
kink crystal^. 

In this symmetry breaking pattern, two rotational and 
one translational symmetries are spontaneously broken. 
We note that there appears no gapless mode associated 
with the broken rotational symmetry. In general, the 
number of NG modes does not coincide with that of bro¬ 
ken global spacetime symmetries [40-42]. 

In the remainder of this section, we present a general 
discussion on the effective low-energy theory of phonons 
(see also [25, 31, 32, 43, 44[). Let us consider a the¬ 
ory with the free energy F[(j)\ = J d^a; and 

assume that {(p) = 4 >q{x) is a static solution minimiz¬ 
ing F[(j)]. In the following, we assume that the free 
energy density T respects rotational and translational 
symmetries. If (j)Q is modulated in one direction and 
is homogeneous in the transverse directions, one can 
write it as 4>o{x) oc f{q ■ x), where q is a vector paral¬ 
lel to the modulated direction and /(•) is a dimension¬ 
less function. Now we consider a translational fluctua¬ 
tion corresponding to the phonon around this solution, 
4'{x) cx f{q ■ X + qu{x)) with q= \q\. Plugging this func¬ 
tion into F[(j)] and expanding in powers of u and Vu, one 
obtains the effective theory for the u field. 


Rotational symmetry of the original free energy implies 
that f{q ■ x) and f{q' ■ x), with a rotated vector q', have 
the same value of ^[ 15 ^]. Writing q' x = q x + (q'~q) x, 
it follows that the fluctuation qu{x) = {q' — q) x does not 
cost any additional free energy. This imposes a constraint 
on the form of the effective theory. Suppose q points 
in the z direction, with no loss of generality. Under a 
rotation by 9 about the y axis, the gradients of u are 
given by 


dzU = —1 -|- cos 9 , dxU = sin 0 , (2) 

hence {1 + dzU)^+ ±u)^ is invariant under the rotation. 

Then the effective theory of u invariant under rotation 
can be constructed as a function of {1+dzu)'^+ {W±u)'^ = 
2dzU -|- (Vu)^ -|- 1 and higher-order derivative terms such 
as 

A total derivative term dzU would be allowed in the ef¬ 
fective theory if we only require the stationary condition 
for F under a variation of (jj with fixed boundary condi¬ 
tions. However, global minimization of F excludes this 
term, which can be seen as follows. Gonsider a fluctua¬ 
tion u = ez over the ground state (j)o{z). It corresponds 
to a dilatation z —)■ (1 -|- e)z. In the presence of a term 
cx AdzU in F, this u will generate an energy shift 

AF = eA + 0{e‘^), (3) 


indicating that there exists a lower-energy direction when 
A ^ 0. In other words, we have A = 0 when the conden¬ 
sate is the ground state. 

This observation imposes a constraint that F should 
depend on 2dzU + (Vrt)^ at least quadratically. Thus 
we conclude that the free energy density for the phonon 
assumes a form 


F = B 


dzU + 


+ c(viw)2 


(4) 


at the leading order of the derivative expansion, with low- 
energy constants B and C. In deriving (4) we exploited 
the fact that the term {dzu){'V']_u) is prohibited if we 
assume parity invariance at low energy (see Appendix 
E 3 for more details). Here we emphasize that the ab¬ 
sence of (V_Lri)^ in F makes the dispersion of u strongly 
anisotropic. 

Such an anisotropic dispersion does not appear when 
the rotational symmetry is explicitly broken, just as in 
QGD under external magnetic fields. This is shown in 
Appendix A. 


^ The free energy has no Lorentz symmetry because it is broken by 
the existence of matter. Therefore, we do not take into account 
the Lorentz symmetry breaking as the breaking pattern. 

^ In general, real kink crystals can be classified by the Frieze group. 
The real kink crystal discussed in this paper has the Z x Z2 sym¬ 
metry, where Z and Z2 represent the glide reflection symmetry 
and the reflection symmetry at a certain vertical line, respec¬ 
tively. 


III. REAL KINK CRYSTAL 

In this section, we study the long-wavelength fluctu¬ 
ations of the real kink crystal chiral condensate within 
the NJL model in the chiral limit. The real kink crys¬ 
tal is known to be energetically favored over a spatially 
homogeneous chiral condensate in a certain range of tem¬ 
perature and chemical potential [15[. We analyze the GL 
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free energy of the NJL model near the Lifshitz point and 
derive the elastic free energy of phonons and pions on the 
kink crystal on the basis of Bloch’s theorem for particles 
moving in a periodic potential. For simplicity we will not 
consider dynamical aspects of those fluctuations and also 
ignore gapped fluctuations. 


A. Phase diagram 

According to Nickel’s work [14], the GL expansion for 
the (3 + l)-dimensional NJL model in the chiral limit, to 
sixth order, reads 

nGL[M(x)] = 02 m 2 + a4 {M^ + (VM)2} 

+ 06 {2M® + 10M2(VM)2 + (AM)2} . 

(5) 

Here the real-valued held Mix) represents the chi¬ 
ral condensate. For simplicity, the pionic condensates 
are suppressed at this stage. The neutral pion 
fluctuation will be later incorporated in Sec. HID, and 
a fully isospin-symmetric form of the GL expansion is 
presented in Appendix B. 

The coefflcients in the GL expansion (5) are given as 
functions of temperature (T) and chemical potential (/r) 
[14] 


02 = 


1 

4G 


04 = 


O 4 


06 = 


12 ’ 


4 = (-l)t41V,iV/T^ [ 

n “'■'sg 


d^p 


( 2 ^)^ [(a;„+i^)2+p2]"/2 


( 6 ) 


where G is the four-fermion coupling in the NJL model 
and Nc and Nf are the numbers of colors and flavors, 
respectively. The momentum integrals for 02 and 04 are 
ultraviolet divergent and must be regularized with some 
regularization scheme. This GL expansion is valid near 
the QGD critical point at which 02 = 04 = 0. We require 
06 > 0 for stability. 

To And the correct ground state of the GL free energy, 
we have to solve the GL equation: 

= (7) 


If we assume a one-dimensional modulation in the z direc¬ 
tion, the GL equation reduces to that of the GN 2 model, 
which has a family of solutions [22-24] given by 

Moiz) = qy^sii{qz;v), ( 8 ) 


where g is a function of iz through the relation 
ly +1 04 , 1 02 


4 , ly -r ± <^4 2 , 

q -|- ——— — q — ■' — 

+ diy + 1 ae ly^ + diy + I 


= 0 . 


(9) 


One can easily check that Mq with (9) satisfies the GL 
equation (7).^ The elliptic parameter iz G [0,1] is not de¬ 
termined by the GL equation alone and must be fixed 
from the requirement of lowest energy per period, ly 
controls the shape of the solution; sn( 2 ;; 0 ) = sin 2 and 
sn( 2 ;; 1) = tanhz. We can identify the solution sn(z; 1) 
as the homogeneous solution because thermodynamically 
tanh 2 : is equivalent to a spatially homogeneous configura¬ 
tion. Let L be the period of Mg and Q the wave number 
of Mg, respectively given by 


4Kfy) 

q 


and 



( 10 ) 


where K(zz) = dt (1 — jz sin^ is the complete 

elliptic integral of the first kind. In what follows, we 
denote the average of a periodic function as 

F=^J^ dzF{z). (11) 

To map out the phase diagram near the tricritical 
point, we need to evaluate ( 6 ) for given T and p. In 
this work, the Matsubara sum in ( 6 ) was done analyti¬ 
cally, whereas the momentum integral was regularized by 
a three-momentum cutoff A. (See Appendix G for addi¬ 
tional comments on the choice of regularization.) In this 
work, we have used the parameter set in [46] : A = 632 
MeV and GA^ = 2.173. The value of jz at a given {T,p) 
was numerically determined through minimization of the 
GL free energy per period 

^Hgl[Mo(z)]. (12) 


Then q immediately follows from v via (9). 

In Fig. 1, we show the resulting phase diagram of 
the GL-expanded NJL model (5). Strictly speaking, the 
GL analysis is only valid in a region with a small or¬ 
der parameter and small spatial variation, so the phase 
structure shown here is only intended to be a qualitative 
guide. In the red region, the kink crystal chiral conden¬ 
sate (0 < iz < 1) is energetically favored, while in the 
green region, the homogeneous chiral condensate (y = 1) 
develops. The three phases (real kink crystal, homoge¬ 
neous and symmetric phases) meet at the Lifshitz point, 
(T,/r) = (81.4,276.4), at which 02 = 0:4 = 0. The real 
kink crystal phase is very narrow near the Lifshitz point 
and is almost invisible in the plot. The phase structure 
from our numerical calculation is consistent with anal¬ 
yses with a non-expanded effective potential [15] except 


® The most general solution to (7) is given by a quasi-periodic func¬ 
tion written by the Riemann theta function with genus g = 3 [45], 
which is considerably more complicated than the elliptic func¬ 
tion, (8). In the following, we assume that (8) is energetically 
favored over higher genus functions, and leave a more compre¬ 
hensive analysis to future work. We thank D. A. Takahashi for 
useful comments on this point. 
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FIG. 1. Phase diagram of the NJL model in the chiral limit 
with 6th order GL expansion. The Lifshitz point is located 
at = (81.4,276.4) [MeV]. 


that the real kink crystal phase in Fig. 1 appears to be 
broader at low temperature. The phase boundary in our 
results are consistent with Nickel’s observation [14] that 
the real kink crystal phase is favored if a 2 > 0 and 



Ct 2 CtQ < 0:4 < 


— 2,^02^ ■ 


(13) 


In Fig. 2, we plot the root-mean-square condensate 





(14) 


and the wave number Q of the real kink crystal at T = 70 
MeV and T = 10 MeV. It is observed that both phase 
transitions (from the kink crystal phase to the homoge¬ 
neous and to the symmetric phase) are of second order: 
At the transition with lower fj., the wave number grad¬ 
ually rises from zero, implying the formation of widely 
separated domain walls. On the other hand, at the transi¬ 
tion with higher n, the amplitude of M vanishes smoothly 
with keeping a nonzero wave number. This behavior is 
consistent with the preceding work with a non-expanded 
potential [15]. 



FIG. 2. The average magnitude of M and its wave number 
Q, at T = 70 MeV (top) and T = 10 MeV (bottom). 


into (5) and expand in u, dropping total derivatives. 
Then 


B. Phonons 


We would like to derive the elastic free energy of 
phonons originating from the spontaneous breaking of 
translation symmetry in the real kink crystal phase. Let 
us substitute 

M{x) = Mo(z -|- u{x)) 

= Moiz) + Mq{z)u{x) + ]^MQ{z)u{xf‘ H- 


nGL[M{x)] 

= nGL[M„(*)] + 

+ hi{z){dzu)(y]_u) 

+ h2iz){d'^u){'\/]_u) + 0(u^) , 
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FIG. 3. The order parameter Mo and the coefficient fnnc- 
tions /i and g\ over one period at (T, /i) = (70, 286.0) [MeV] 
(top) and = (70,286.5) [MeV] (bottom). All functions 
are measured in units of A. The period L is 7.71 fm (top) 
and 6.40 fm (bottom), respectively. 


where Vx = {dx,dy) is a transverse derivative and fi, 
/ 2 ) ffi) 92 , hi, /i 2 are defined as 

fiiz) = 2ia4 + 10aeM^){M^f 

+ 4a6{(M")"-2M'M'"} , 

/ 2 (z) = 2 a 6 (M') 2 , 

gi{z) = 2(^4 + 10a6M2)(M')2 - dagM'M'" , 

52(2) = 2a6(M')2, 

hiiz) = dagM'M" , h2{z) = 2a6(M')2 . 

Figure 3 displays the order parameter and some of the 
coefficient functions at (T,fi) = (70,286.0) [MeV] and 
(T, fi) = (70, 286.5) [MeV]. Obviously the functions share 
the same period with the order parameter. Note that 
(Vxw)^ is present in (15) even though the general theory 
in Sec. II suggests that this term should be absent in the 
effective action of low-energy modes. This is not a con¬ 


tradiction: since the expansion (15) includes functions 
/i, gi, etc., that vary over the microscopic length scale 
L, we cannot immediately infer the dispersion of modes 
with wavelength much longer than L from there. 

As emphasized in Sec. II the free energy of a modulated 
condensate must be invariant under a spatial rotation. 
This can be checked for (15) as follows. An infinitesimal 
rotation about y axis is equivalent to a displacement field 
u{x) = — sin 0x-|-(l —cos 0 ) 2 : with |0| ^ 1. Plugging this 
into (15) and expanding in 9, we find the leading term 
to be hg^(^z)9'^. Then, for the rotational symmetry to be 
preserved, the average of gi must vanish: 


j gi{z) = 0 . 


(17) 


A direct proof of this equality based on the minimization 
of energy is given in Appendix D for completeness. We 
remark that (17) is not automatically ensured by the GL 
equation (7) [or equivalently (9)] alone—we must min¬ 
imize the energy per period, to have (17) satisfied. It 
will be shown below that the property (17) is instrumen¬ 
tal in making the dispersion of phonons anisotropic, in 
accordance with the general argument in Sec. II. 

To evaluate the low-energy phonon fluctuation on the 
real kink crystal, let us consider the eigenvalue equation 


Eu = 


(5F[u] 

Su 


(18) 


with 


F[it] = J 

51W 




. 2 
(Vxu)" 


92{z) 


{Vluf 


2 ^ ' 2 
hi{z){d^u){V\u) + h 2 {z){dlu){V\u) 


(19) 


The eigenvalue equation in the explicit form reads 

HuU = Eu (20) 


with 

Hu = -d.ifid.) + dlihdl) - giVi + 52Vl 

-{d.hi)wl + wl{h2,dl]^, 

where denotes the anti-commutation relation and 
{h2,dl) _i_ u = h 2 d‘^u -I- dl{h 2 u). As the operator Hu 
acting on u is real and Hermitian, the eigenvalue E is 
real, and we can elevate u to a complex-valued function 
without changing the eigenvalues. We note that E it¬ 
self does not give the dispersion relation of the phonon, 
but is rather related to the phonon susceptibility. The 
dispersion relation can in principle be obtained from the 
time evolution equation, but it is complicated because 
of the mixing with hydrodynamic modes (see, e.g., [32]) 
and will not be discussed further in this paper. 

Since all the coefficient functions are periodic functions 
sharing the same period, it follows from Bloch’s theorem 
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FIG. 4. Eigenvalues of Hu in (21) at = (70,286.0) 

[MeV], for = 0 (top) and for = 0 (bottom). The 
domain —0.5 < k^/Q < 0.5 is the first Brillouin zone, kz 
and E are normalized by Q, and k± is normalized by the UV 
cutoff scale A. 


FIG. 5. B and C at T = 70 MeV (top) and T = 10 MeV 
(bottom). B and C are normalized by the UV cutoff param¬ 
eter A with appropriate dimensions. The dotted vertical lines 
mark the boundaries of the modulated phase. 


that u can be decomposed into a plane wave and a peri¬ 
odic function, 

u(a:) (22) 

where fc_L is the momentum in transverse directions, 
is the so-called crystal momentum, and (j){z) is a periodic 
function, viz. (j){z + L) = (j){z). Substituting (22) into 
(20) yields an eigenvalue equation for cj), which we have 
solved numerically by way of a Fourier decomposition 

= J2n=-n„,uz "-“ax = 20. To 

see convergence, we have increased rimax up to 30 and 
confirmed that the results are unchanged. 

In Fig. 4, we show the eigenvalue E numerically com¬ 
puted for varying kz and k±. A marked difference from 
the eigenvalue of particles in a free space is that there 
are infinitely many levels for given momenta, in analogy 
to electrons in metals which develop a band structure. 

It is the lowest eigenvalue Eq (red curves in Fig. 4) that 


pertains to the free energy of long-wavelength phonons. 
By adopting a variational approach, one can rigorously 
show that Eq behaves for ^ A:_l ~ 0 as 

Eq ~ Bkl + Cki , (23) 

where B and C are functions of T and yL. The absence 
of the 0{k^) term in (23) is guaranteed by the property 
(17). The proof of (23) is somewhat technical and is 
relegated to Appendix E. Equation (23) shows that the 
elastic free energy of low-energy phonons becomes 

Ki = lJ d'a: [Bidzuy + C{Wluf] . (24) 

One may suspect that the coefficient B would be given 
by § fi- However, as shown in Appendix E, this naive 
guess is incorrect; the coupling between (pQ and (pn^o is 
not negligible even in the perturbation series in kz- 
To extract the values of B and C from eigenvalues, we 
have numerically fitted the curve of Eq with trial func- 
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tions Eq = Bk^ and Eq = Ck^, respectively, with B and 
C treated as fitting parameters. Figure 5 shows B and 
C at T = 70 MeV and T = 10 MeV obtained this way. 
Notably, B and C in Fig. 5 are positive throughout the 
real kink crystal phase, which proves local stability of this 
condensate in agreement with numerical results in [17]. 
Because the phonon mode exists only in the real kink 
crystal phase, it is natural that both coefficients tend to 
zero at the phase boundaries, although the eigenvalues E 
were too small near the left boundary for T = 70 MeV 
to perform a reliable fitting. 

We mention that the spectrum of gapless excitations 
over the same background (8) has also been worked out 
in [45] with entirely different methods. However a direct 
comparison is difficult owing to the fact that the model in 
[45] is non-relativistic and in one space dimension, while 
our model is relativistic and in three space dimensions. 


C. IR divergence and quasi-long-range order 


Next, we wish to evaluate the impact of thermal fluc¬ 
tuations of phonons on the stability of the real kink 
crystal. Taking the Fourier decomposition M{x) = 
M-n /a/l, treating u in the Gaussian ap¬ 

proximation and ignoring the pion fluctuation, we obtain 




1 

7l 






(25) 


with 


^) = 


27r 

(27r)" 


dfc_L k± / dfc 


T 


^Bkl + Cki 


T , i± 
;log 


(26) 


AttVW ^/cJb ' 


where we have only incorporated the lowest Matsubara 
mode since it is dominant in the infrared. The momen¬ 
tum integral is IR divergent and is regularized by a cutoff 
t'_L. One can regard £_l as the transverse diameter of the 
quark matter in a compact star. In the thermodynamic 
limit ^_L —>■ oo, the condensate (25) drops to zero with 
negative powers of implying that the one-dimensional 
modulation is wiped out by thermal fluctuations at any 
low T > 0, a phenomenon known as the Landau-Peierls 
instability. In fact, if the average amplitude of displace¬ 
ment fluctuation exceeds the interval of layers, it does not 
make sense to speak of a spatial long-range order. We 
emphasize that this instability persists even at nonzero 
quark masses, since it originates from phonons that re¬ 
main elastic regardless of the quark masses. 

In the thermodynamic limit, the system instead ex¬ 
hibits a quasi-long-range order. Expanding M in Fourier 
series and ignoring non-Gaussian effects and pion fluctu¬ 
ations, we find that the correlation function of the order 



p [MeV] 



FIG. 6. The exponent characterizing the algebraic decay 
of the order parameter correlation function in the kink crystal 
phase at T = 70 MeV (top) and T = 10 MeV (bottom). The 
dotted vertical lines mark the boundaries of the modulated 
phase. 


parameter behaves for |a:| L as 
{M{x)M{0)) 

= E 


_ -Mn^m ^inQz {nu{x) TOu(O))]) 

Lj 

n,m 


= E 


MeMel exp 

L 


— ^ (^{nu{x) -\- mit(O))^) 


^ M„Mm ^inQ. 


= ^ ^ g.gjj^p 


Q'^ f (rk ni^ + + 2mne^'^'^ 


-T 


2 J (27r)3 -h Cki 


= exp 


-n^Q^T, 


d-^k 1 - e*' 


E 

n> 1 


\Mn\'^ / 2cos(nQz) x jzj {x± = 0) . . 

nA ■ ( ) 


(27r)3 Bkl E Ck\ 

[x^ = { 

(z = 0) 


In the intermediate step, we have dropped terms with 





































n ^ —m since their momentum integrals are infrared 
divergent. The exponent r]c above is defined by 


^ Q^T 
~ SttVW 


(28) 


which was originally introduced by Caille [47]. Such an 
algebraic decay of the order parameter correlation has 
been known to appear in the smectic-A phase of liquid 
crystals [25, 31-33], in the FFLO phases of fermionic su- 
perfiuids [34, 37, 38] and in the modulated pion conden¬ 
sation in nuclear matter [48] . Such a slow decay of the 
order parameter suggests that it would be hard in prac¬ 
tice to distinguish it from a true long-range order. 

In Fig. 6, we show the critical exponent rjc at T = 70 
MeV and T = 10 MeV. Within numerical precision, rjc ap¬ 
pears to go to infinity at the ends of the modulated phase. 
This is physically acceptable because an exponential de¬ 
cay of the connected two-point function is expected both 
in the symmetric phase and in the homogeneous broken 
phase. 

Although {M(x)) = 0, a closer look at steps leading to 
(27) shows that {M{x)‘^) yf 0 in this phase. The real kink 
crystal phase is thus characterized by a homogeneous 
higher-order condensate consisting of four quarks, which 
is qualitatively distinct from the naive mean-field phase 
with {M{x)) 0. These phases may be distinguished 

by (Z 2 ) R X (I-2)l symmetry that flips signs of quarks of 
each chirality independently. Such a novel higher-order 
condensate is discussed for the FFLO phase of nonrela- 
tivistic fermions in [34, 37, 38] . 

The large fluctuation of phonons can be suppressed 
by several factors: (I) strictly zero temperature, (II) 
higher-dimensional modulations, (III) coupling to an ex¬ 
ternal vector field (e.g., a magnetic field) and (IV) a 
finite volume. In case II, the phonon fluctuations no 
longer cause infrared singularity because the number of 
spatial directions with a quadratic dispersion decreases 
[25]. However, it is still an open problem whether higher¬ 
dimensional modulations can be energetically favored 
over a quasi-long-range-ordered phase. Case III is due 
to the fact that, as laid out in Appendix A, the explicit 
breaking of rotational invariance leads to a non-vanishing 
(V_lu)^ term. This leads to an interesting observation 
that inhomogeneous condensates in QCD under magnetic 
fields [49-53] could be stable against fluctuations. 

Next, we shall analyze the finite-volume effect of case 
IV in detail. If the system size is finite, (u^) remains 
finite because the IR divergence is cut off. This is true 
for putative quark matter inside neutron stars. Roughly 
speaking, if (u^) < the one-dimensional structure is 
expected to remain, while it is likely to be wiped out 
when (u^) > L^. Using (26), we define the crossover 
length as a scale below which the one-dimensional 
modulation persists. 


a = \/C/B , (29) 


which grows rapidly as the temperature decreases. The 



[MeV] 



FIG. 7. The crossover scale and the period L of the 
kink crystal at T = 70 MeV (top) and T = 10 MeV (bot¬ 
tom). The dotted vertical lines mark the boundaries of the 
modulated phase. Although L = 2ti/Q diverges at the left 
boundary of the modulated phase (recall Fig. 2), this is not 
visible in these figures due to limited numerical precision. 


periodic structure of the condensate will survive if the 
size i of the quark matter is in the window L ^ 

Figure 7 shows the crossover scale ^_l together with 
the period L of the real kink crystal at T = 70 MeV 
and r = 10 MeV. It is observed that L is of order 1 fm 
throughout the modulated phase, while ^_l has a strong T 
and /i-dependence. At high T, soft phonons are so easily 
excitable that the crossover length is comparable with 
the period. Thus, there is no vestige of a one-dimensional 
crystal at high temperature. On the other hand, at low 
T (< 10 MeV), the situation is different. As shown in 
the lower panel of Fig. 7, there is a region in which is 
macroscopic, i.e., of order 1 km (= 10^® fm), which would 
be presumably large enough to accommodate a quark 
core of neutron stars. Although the GL expansion is not 
quantitatively reliable at such a low temperature, our 
conclusion that the real kink crystal condensate persists 
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only at very low temperatures should be qualitatively 
correct. 


D. Pions 

Next we proceed to the analysis of gapless pions that 
stem from the spontaneous breaking of SU(2)fl; x SU(2 )l 
to SU(2)y. Unlike phonons that only emerge in the mod¬ 
ulated phase, pions show up both in the homogeneous 
broken phase and in the real kink crystal phase, and 
hence it is of great interest to investigate the nature of 
pion fluctuations across the transition between these two 
phases. Since the vectorial isospin symmetry is intact 
in both phases we will only take account of the neutral 
pion fluctuation for simplicity, and defer a fuller analysis 
to Appendix B. As in [14, 15], let us make M complex 
as M{x) = —2G[S'(a:) -I- iP 3 {x)], with S{x) = {'il!tp{x)) 
and Psix) = {tpijsT^iIjix)). It has been shown by Nickel 
[14, 15[ that, assuming a one-dimensional modulation, 
the (3 -I- l)-dimensional GL Lagrangian for M can be 
obtained from that of the chiral Gross-Neveu model [22- 
24[, leading to the result 

UGL[M(a:)] = a2|M|2 + a4 {|M|4 + |VM|2} 

-ka6{2|M|®-h8|M|V^I^ (30) 

-k 2 Re -H IAM|2} . 

Now we shall follow the same route as for phonons. In 
the following, tto will be denoted as tt in order not to 
clutter notation. Substituting 

M(a;) =Mo(z)e*"(“) (31) 

into (5) and expanding in tt, we obtain 

[Af(a:)] 

= Ugl[Mo(z)] + 

+ ^i^(Vj_^)2 + ^^(Vi7r)2 (32) 

+ hi^(z){d^Tr){W']_Tr) 

+ /i2^(z)(a^7r)(V5_7r) -h 0{Tr^) , 

where /i^, / 2 ^, 5 i^, 527 r, /i 27 r are defined as 

h^{z) = 2 iai + 6 aeM^){Mo)^ 

+ 4a6{2(M')^-MoM"} , 

/27r(^) = 2a6(Mo)^ , 

= 2(^4 + )('^o)^ ~ ‘^CXqMqMq , 

92 - 17 ( 2 ) = 2 ae(Mo)‘^ , 

hin(z) = 4:aQMoMQ , h 2 -n-iz) = 2 aQ(Mo)'^ . 

We note that, unlike the phonon, 

^ 51. ^ 0 . (34) 


The eigenvalue of pions may be derived from the 
eigenvalue equation 


Ett-k = 


SF[ 


TT 


with 


Stt 

n 2 , / 2 .(z), 


(35) 


F[tt] ^ Id^X {^^id.TT)^ + ^( 5 »^ 
+ ^(V^7r)2 + ^(Vi ^)2 


(36) 


+ hi^{z){dzTT)(\7']_TT) + /i2.(2:)(5f7r)(vi7r)|. 

In a more explicit form, it reads 

Ht^tt = Et^tt (37) 


with 

77. = -d.ifi.d,) + dl(f2.dl) - 5i^Vi + 52 .Vi 

-(a./ri.)Vi+Vi{/i2.,9,^}+. ^ ^ 

Since the structure of (37) is identical to the case of 
phonons ( 20 ), one can apply the same techniques to solve 
it. On the basis of Bloch’s theorem, tt may be decom¬ 
posed as 

TT(x) = e^'^^-^^F^’‘^ (l)(z), (39) 


for a crystal momentum kz and transverse momenta fcj^. 
4>(z) is a periodic function with period L. Substituting 
this decomposition, we arrive at an eigenvalue equation 
for (j), which can be solved with the same numerical meth¬ 
ods as for phonons. 

In Fig. 8 , we show the eigenvalues for varying kz 
and A:j_. To analyze the lowest eigenvalue E^^^q (red curves 
in Fig. 8 ), we have used a variational method along the 
lines of Appendix E, which shows that the leading be¬ 
havior of Ett^ near = fc_L = 0 is given by 

E^^o ~ F^^kl + Fikl , (40) 

where 


Fl = 



(41) 


and is a positive function that has to be computed 
numerically. The elastic free energy of low-energy pions 
is therefore 


EJ = / d^a 


F|f(5,^)2+Fl(V^7r)' 


(42) 


Since Fj ^ 0, pions on the real kink crystal have a lin¬ 
ear dispersion in all directions, in contrast to phonons. 
Thus the thermal fluctuations of pions neither cause any 
infrared divergence nor destroy the long-range order. 

FJ and F|^ at F = 70 MeV are shown in Fig. 9. In 
the homogeneous broken phase, the system is isotropic 
and F]_ = F^. In the kink crystal phase, the system is 
anisotropic and the coefficients no longer coincide. In our 
calculation, in the kink crystal phase, F| > F^ holds at 
any temperature and chemical potential. 
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FIG. 8. Eigenvalues of H-^ in (38) at (T, /j.) = (70, 286.0) 
[MeV], for = 0 (top) and for = 0 (bottom). The 
domain —0.5 < k^/Q < 0.5 is the first Brillouin zone, kz and 
k± are normalized by Q and A, respectively. 


IV. CONCLUSION 

In this paper, we presented a first systematic study of 
low-energy fluctuations in the real kink crystal phase of 
the NJL model. The elastic free energy for the phonon 
was shown to be (9zu)^ in the longitudinal direction and 
{d^uY in the transverse directions with respect to the 
modulation of the condensate, which had an important 
consequence of vanishing order parameter, exhibiting a 
quasi-long-range order. We argued that, since the cor¬ 
relation function decays only algebraically, the real kink 
crystal may be sustained in compact star cores if the tem¬ 
perature is sufficiently low. Our analysis, which should 
be reliable at least near the critical point, suggests that 
fluctuation effects that are missed in the mean-field treat¬ 
ment change the nature of inhomogeneous chiral conden¬ 
sation qualitatively. 

This work can be extended in various directions. The 
NJL model used here can be extended with the inclusion 


FIG. 9. Fl and at T = 70 MeV (top) and T = 10 MeV 

(bottom). 


of vector interaction [16], or one can use the quark-meson 
model [15, 19[; in either case, the phase structure would 
be quantitatively modified. For a fuller understanding of 
the phonon fluctuation, we need to consider anharmonic 
effects due to the O(u^) term in (4), which is known for 
the case of smectic liquid crystals to modify the scaling 
behavior (27) [54, 55[. It is also intriguing to study the 
interaction between quarks, phonons and pions. As for 
the fate of the Lifshitz critical point in Fig. 1, we point 
out that it must be eliminated from the QCD phase di¬ 
agram once fluctuations are fully incorporated, because 
the lower critical dimension of the isotropic Lifshitz crit¬ 
ical behavior with continuous symmetry is 4 [56[. This is 
true in the chiral limit, while for nonzero quark masses, 
the continuous symmetry is broken, and a more careful 
study is needed. At any rate, the mean-field picture can 
fail even qualitatively near the Lifshitz point, and it is 
highly desirable to develop a method for analyzing fluc¬ 
tuations that goes beyond the GL expansion. This is a 
hard problem at this stage but should be seriously con¬ 
sidered in future research. 
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Finally we note that for phenomenological applica¬ 
tions to the physics of compact stars it would be im¬ 
portant to incorporate effects of nonzero quark masses, 
isospin chemical potential, electromagnetic interactions 
and color superconductivity, which deserves further in¬ 
vestigation. 

Note added 

While this work was being completed, we learned of an 
independent work [57] where the low-energy fluctuations 
over a Fulde-Ferrell type inhomogeneous chiral conden¬ 
sate were discussed. That condensate breaks transla¬ 
tional and internal symmetries in a different way from 
the real kink crystal considered in our work, and hence 
produces a different number of Nambu-Goldstone modes. 
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Appendix A: Phonon effective theory coupled to a 
vector field 

In this Appendix we present a simple argument based 
on the method of [44] showing that a coupling to an ex¬ 
ternal vector held modifies the dispersion of phonons in 
a qualitative manner. 

Let us start with a theory with no vector held. Sup¬ 
pose {(p) = (/'o('Z) is a modulated static solution that 
minimizes the rotationally symmetric free energy F']^] = 
/ d^a: T{4i,d4>). Now we consider a translational fluctu¬ 
ation corresponding to the phonon around this solution, 
(p^x) = cj)Q{z + u{x)). Plugging this into Flp] and ex¬ 
panding in powers of u and Vu, one obtains the effective 
theory for the u field. Since u always appears in the form, 
z + u, the effective free energy is constructed from 

scalar functions of z + u, their derivatives, and 6ij. 

(Al) 

More concretely, the ingredients with one and two deriva¬ 
tives are 


di{z + u) = S- + diU , (A2) 

didj{z + u) = didjU , (A3) 

respectively. These ingredients are compatible with the 
original symmetry as they should be. For example, 

(S^ + d^uf = 1 + 2d,u + (Vu)2 (A4) 


combine these ingredients and derive the leading phonon 
free energy 


F = A 




B 


n 2 


dzU + 


(AS) 


Just as we discussed in Sec. II, the total derivative term, 
dzU, linear in u is prohibited by the minimum-energy 
condition, so that we have A = Q and F reduces to (4) 
after taking into account higher derivative terms. The 
reader is referred to [44] for more details of this method. 

We shall then proceed to a discussion on the elastic free 
energy in the presence of an external vector held u*. Such 
a situation is pertinent to modulated chiral condensates 
in QGD with an external magnetic held. Now one can use 
u® in addition to the previous ingredients (Al), so that 
the following interaction can appear in the free energy 
for example: 


Vi (Jf -I- diu) =v''+ v’‘diU . 

The simplest modification of (A5) will then be 


F = A 


dzU + ^(Vu) 


+ B 


dzU F 


(A6) 


(A7) 


-v^ - FdiU. 


We notice here that the total derivative terms in the lin¬ 
ear order arise from both of the A and u® couplings. If 
the vector Vi is an external held, i.e., a non-dynamical 
held, and takes a value 

u® = ASl , (A8) 

the linear order terms cancel out: 

jz= + B(dzuf+ 0(u^). (A9) 

Therefore, when an external held u® explicitly breaks the 
rotation symmetry, the (V_lm)^ term need not vanish, 
which can be realized without spoiling the minimum- 
energy condition of the condensate. The dispersion 
of phonons becomes linear in all directions, implying 
that the severe infrared divergence at finite temperature 
(cf. Sec. IIIG) is ameliorated. 

It would be important to note that the non-vanishing 
(V_lu)^ term in the free energy does not arise if the vec¬ 
tor is dynamical (i.e., not external) and its condensation 
is aligned in the z-direction. In such a case, we have 
gapped fluctuations associated with the spontaneously 
broken rotational symmetries, just as in the smectic-A 
phase of liquid crystals. The free energy after integrating 
out those gapped fluctuations turns out to be the same 
as the one without dynamical vector fields, (4), which 
results in a strongly anisotropic dispersion of phonons. 


Appendix B: GL expansion with pions 


automatically reproduces the rotationally invariant com- In this Appendix, we incorporate pionic modes into 
bination derived in Sec. 11. It is now straightforward to the effective theory (30) so that it becomes manifestly 
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invariant under SU(2)j{ x SU(2)i. It is convenient to Thus the desired GL function with manifest chiral sym- 
work with the matrix field metry is given by 


S(a:) = -2G[S{x)l + iPa{x)T^], (Bl) 

where {r“} are the Pauli matrices. Under Ur G SU(2)i^ 
and Ul € SU(2)j;,, it transforms as S —)■ Ur'SUI^. Then 
the most generic GL function invariant under SU(2)^ x 
SU(2)i is given, up to 6 th order in fields and derivatives, 

by 

Ugl(S) = y trptq + ^tr[a,Eta^s] + ^ (trptS])' 
+ ^(tr[EtE])%|itr[AStAE] 

+ /3i trpl'E] tr[5,E^9,S] 

+ / 32 {tr[( 5 ,E)E^( 9 ,E)S^]+h.c.} 

+ P 3 (tr[E^ 5 iE])%/? 4 tr[E^E] tr[E^AE], 

(B2) 

where the coefficients /3i, f 32 , fis and are yet to be 
determined. However, the last four terms of (B2) are 
not independent, owing to the relations 

tr[(9,E)Et(5,E)Et] 

2 1 . . (B3) 

= (tr[Et5,E]) - - trp^E] tr[a,Eta,E], 

and 

tr[El'E] tr[E^AE] 

= d, (tr[EtE] tr[Eta,E]) - 2 (trptS.E])" (B4) 
-tr[E^E] tr[5iE^5iE] . 


Therefore one can let P 2 = Pa = ^ without loss of gener¬ 
ality (up to total derivatives). Since Ugl(E) must reduce 
to (30) when Pi = P 2 = 0, we require 


f^GL(E) 

= ^ tr[EtE] + ^ tr[5,Et5,E] + ^ (trptE])' 

+ ^ (tr[EtE])^ + ^ tr[AEt AE] 

+ ^ tr[EtE] tr[a,Eta,E] + ag (trpt^.E])" . 

Appendix C: Regularization of the thermodynamic 
potential 

In the main body of this paper, we have employed the 
momentum cutoff to regularize UV divergences in the 
NJL model. Although the momentum cutoff is subtle in 
the presence of a modulated condensate as is claimed in 
[14, 20], our work is concerned with the GL expansion 
around a homogeneous vacuum, and so this problem is 
unlikely to obstruct our present analysis. We would also 
like to add that the phase structure including a modu¬ 
lated phase has been found to be quite robust against 
varying regularizations [58]. That being said, however, it 
would be desirable for theoretical completeness to have 
an alternative derivation of regularized GL coefficients 
that is free from the above subtlety. In this Appendix we 
shall demonstrate one of such regularization methods. 

Let us start from the mean-field thermodynamic po¬ 
tential of the NJL model in Euclidean spacetime, 

= -N^Nf log det [0 - ^74 -f M(a;)] 

+ /dS ^ 

The functional determinant suffers from UV diver¬ 
gences, and one has to specify a regularization scheme. 
The conventional three-momentum cutoff is not useful 
for a generic inhomogeneous mean field, whereas the 
Schwinger proper-time regularization [59] is tricky at 
nonzero chemical potential [60, 61] . Now we introduce 
a new regularization scheme free from these problems, 
based on the formula 


4 / 3 i|MnVMp-K 2/?3 {iMplVM^-k Re [(VM)^!^*^] } . 

3 

.-. Pi = -06 and P 3 = aQ. (B5) 


/■°° f 1 11 

Jo l(i« + A)2 (^ + 5)2/ 

(G2) 

See [62, 63] for the relation of this approach to the naive 
proper-time regularization. The methodology of deriva¬ 
tive expansion is well known [64], and we shall be brief 
here. 
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log = I - A + (fM) + M^] - log [-(94 - - A] } 






X J p J 0 


dwwe 


dww 


1 


1 


X J p J 0 


[w - (94 - /i)2 - A + (^M) + M2]2 \w - (94 - ^lY - A]2 
1 1 


[w + (p 4 + — (V + zp)2 + (^M) + M2]2 [ix; + (p4 + i^)2 -(- p2j5 


dw 


2 u> 


X */ p «/ 0 


( (ti; + P2)3 


trO - 


Stw 


(w + P2)4 


trO" 


4w 


(w + P2)5 


trO^ - 


Sw 


(w + P2)6 


trO"* 


where the “tr” denotes a trace over spinor indices, A is 

aUV cutoS, I 

—2ip • V — A + (^M) + and = (p 4 + + p^. 

For the purpose of obtaining the GL expansion up to 6 th 
order, it suffices to expand only up to O'^. After a bit of 
algebra involving a trick piPj —>■ 5ijP^/i and integration 
by parts, we obtain 


tuned that Mq(z) attains the minimum of the free energy 
per period. To show this, we consider a scaled configu¬ 
ration 


M{X,z) = Mo{Xz), (Dl) 


tr 6 = 4M^ , 
trO^ = 4[(VM)2 -hM^], 
tr63 = 4[M®-k7M2(VM)2-K(AM)2], 

- 1 

trO^ = —p2[4M2(VM)2 {AMf ], 

where we have omitted total derivatives, terms that are 
higher order in the GL expansion, and terms odd in 
p. Using a formula ^ Ip 

match the expanded determinant with (5) and extract 
the regularized GL coefficients as 

a 2 

« = 35“ 

a, = -iN,N,IJ^ (C4c) 

In the limit A —>■ oo, we formally recover ( 6 ). This deriva¬ 
tion, in which the functional determinant is regularized 
directly, satisfies the requirement [14, 20] that surface 
terms in momentum integrals strictly vanish. However, 
the three-momentum cutoff is practically more useful and 
is used throughout the main part of this paper. 


Appendix D: Proof of 


51 = 0 


which has a period L/X with L defined in (10). It follows 
from the minimum-energy requirement for Mo{z) that 
the GL free energy per period of M{X,z) must have an 
extremum at A = 1, namely 


r ^ 
Ah 
A—ll OX 


1 

p/A 


/.i/A 

/ dz 

Jo 


Ugl [Af(A, z)] 


= 0 . 


(D2) 


Note that this equality is trivial if iXQi^[M] did not de¬ 
pend on the derivatives of M, for one can simply get rid 
of A from inside of [... ] via a change of variable z —>■ z/A. 
However, this is not possible for U[M] in (5) which does 
depend on the derivatives of M. 

An explicit calculation yields 


A 


1 

lJx 


.L/A 
/ dz 

Jo 


UGL[Af(A, z)] 


= 0 + @ + 0, (D3) 


In this Appendix, we prove (17), i.e., that 5 i(z) must 

vanish in the average sense when iz and q in ( 8 ) are so where [with M' = dzM{X,z) and M” = d‘^M{X, z), and 
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omitting the subscript “GL” for brevity] 


0=1 


^L/A 


dz n[M{X,z)]. 


g = -jn[M{X,L/X)] = -jn[Mo{L)], 
,_, 1 a 


(D4) 

(D5) 

(D6) 


L/X 

I 


f dn dn ^ 


.L/A 

' dz 
0 

dn 

dM' 


{ dn dn 

^ dM' 

M' + 2-^ M” 


n2 

dM" ^ 

dn 


dM 

'W 


(D7) 


dM' 


■M" 


dM" 


1 r dO 


dz 


dn 


dn 


z -^- M' + 2- 

dz dM' dM' 


M" 


1 1 

= -n[Mo{L)] — / dz n[M] 

^ L Jo 


2^ /-L/A 

I L 


i9r2 „ dV, .. 
:M' + 2t7W7tM" 


dM' 


dM" 


(D8) 

(D9) 

(DIO) 


In the step from (D7) to (D8), we used the relation 

^ = ^M'. Now, plugging rfl, ^ and into (D3) 
dX X , , '—^ ^ ^^ 

and taking the limit A —)■ 1, we obtain 


,w „ dfl 
:M' + 


dM' 


dM" 


= 0 . 


(Dll) 


M=Mo 


long-wavelength phonons is given, at leading order, by 
(23). First of all, we note that, among infinitely many 
energy levels that follow from (23), it is only the lowest 
level Eq that matters for the low-energy phonons. Then 
it is easily seen that Eq for given k = {k±, kz) is given by 
the formula 

= .SiS) fTp / 

+ i9i+dzhi)\y^uf+g2\Vluf (El) 

-I- h2{d'^u){'V]_u) + /i2(V^?l)(9fM)|, 


where il(fc) stands for the set of smooth functions of the 
form 


u(a;) <()(z), 

(l){z + L) = (j){z), kz,k±GR. 


(E2) 


1. Eigenvalue with kj_ 7 ^ 0 and kz = 0 

Now we consider Eq for nonzero transverse momenta. 
Substituting (E2) with kz = 0 into (El), we get 


Eo{k±) = min ■ 


1^1^ 


/l|</>'P + /2k"P 


+ (^±51 + ^752)1X1^ ~ k ‘]_ h 2 {( l )'' 4 ' + 4 ‘ 4 ''') ) 


(E3) 


Recalling that flcLiAF] is given by (5), this translates 
into 

^ [2(a4 + 10a6M7)(M')2 + daeiM^f] = 0 , (D12) 

which reduces to ^ = 0 via integration by parts. This 

completes the proof. We stress that the stability of the 
condensate under dilatation has played an essential role 
here. 


where the primes denote derivatives by z and we defined 
gi = gi + h'l. When k± = 0, the minimum is trivially 
Eq = 0, which is achieved by an arbitrary constant solu¬ 
tion: (j){z) = (j)o d- Now for sufficiently small k±, we 
can expand J) that corresponds to the minimum of (E3) 
in perturbative series of k± as 

4>{z) = (j)Q[l + kj_4>i{z) + k\(j) 2 {z) -I- ...] . (E4) 

Plugging this into (E3) yields 


Appendix E: Variational analysis for the lowest 0^o(*_l) = min {j52[(t)i]k\ + I53[4>i]k± + l3A[(t>i]k]_ + 0{k\)] 

eigenvalue spectrum of phonons 

= /32*fc7 0 -P disfkJ^ -p , (E5) 

In this Appendix we apply a variational technique to 
the eigenvalue problem (20) and show that the energy of with 



























15 


I32[Pi] = j 

f {/l '/’iP + f2\PlY + 51 } > 

(E6) 

PalPi] = j 

f {/l(^2<i^’'l + <I>W2) + f2{p2p'i + P'(^2) + + Y>l) - k2{P'{_ + 0l)} - P2[P] j {Pl + ^l) , 

(E7) 

/34[<)>i] = J 

f |/l( </’ 2 p + 4>'i4’3 + + / 2 ( <(> 2 ^ + + '/' 3 V 1 O + 5l( '/’iP + 4’2+ P 2 ) + 92 



- h 2 {p'iPl +Pip'i + P 2 + ?!> 2 )} - ^2[0]^ {\PlY +p 2 +P 2 ) - PsiP] j>{pl +Pl) ■ 

(E8) 


Since EQ{k±) is the minimal eigenvalue, the leading co¬ 
efficient /? 2 * must be at the minimum as a functional of 
(j), i.e., 


P 2 * = min^ {/iP + /2 1*/*! P + 91 } ■ (E9) 

As /i and /2 are positive functions, the minimum trivially 
occurs when =0. Thus, 






(ElO) 


where the last equality is proved in Appendix D. 

Next substituting the solution = 0 into we 

readily find that the coefficient of k\ vanishes: 


, 03 * = min 

0iGC 


+ (pi)} - h* 


(</>i + </>i) 


= 0 . 

(Ell) 


Finally we come to the coefficient of k\_. Since it is the 
leading nonzero term in EQ{k±), it must be minimized 
as a functional of (p. Substituting <p'i = 0 into fii\(pi] and 
using /32* = /Sa* = 0, we obtain 




= ^52 + {^fi\(p'2? + f2\4>2? 

+ gi{.4>2 + <p2) ~ ^ 2(02 + ^ 2 )^ ■ 


(E12) 


Although it is not analytically tractable, one can in prin¬ 
ciple determine p 2 that minimizes the integral by solving 
the Euler-Lagrange equation 

- {M' 2 )' + (/2</>2)" + 51 - K = 0 (E13) 


with periodic boundary conditions. 

Summarizing above, the eigenvalue Eq for transverse 
momenta is given by 


(E14) 


2. Eigenvalue spectrum with ^ 0 and k± = 0 

Next we consider Eq for parallel direction, 

Eo{k,) 

. f {fi\ikz<p + + f 2 \klp - 2ikz<p' - (p"\^] 

= -mF-■ 

(E15) 

When kz = 0, the minimum is trivially Eq = 0, corre¬ 
sponding to a constant solution p{z) = (po 0. Now, for 
sufficiently small kz, one can expand p that achieves the 
minimum of (E15) in a perturbative series of kz as 

p{z) = Pq[ 1 + kzXi{z) + klx 2 {z) + ...] ■ (E16) 

Substituting this expansion, we find 

EQ{kz) = min {x 2 [xi\kl + 0{kl)] , (E17) 

4> 

with 

l 2 [pi] = ^ {/i|l - ix'i? + / 2 |x'iT} ■ (E18) 

Since Eo^kz) is the minimal eigenvalue, xi must be cho¬ 
sen so as to minimize 72 [xi]- Thus, xi satisfies the Euler- 
Lagrange equation 

= -*/( - (/ix'i)' + if2XiY' = 0 , (E19) 

Wi 

with periodic boundary conditions. Using the solution 
Xi*, the eigenvalue at small kz is finally obtained as 

Eoikz) = X2*kl + 0{kl), (E20) 


with 72 * = 72 [Xi*]- 


3. Absence of kzk\ term in Eo 

We have shown the leading behavior of Eq for k± ^ 0, 
kz = 0, and for fc_L 0, = 0 in the previous subsec¬ 

tions. Here we show that Eq does not have a kzk"^ term, 
when both kz and k± are nonzero. 


Eo(fc±) =/34*fcl+0(fci). 
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Since Hu in (21) is invariant under z —)■ —z, the eigen¬ 
values with the momentum kz and —kz degenerate, which 
implies the absence of the kzk^ term. To see this explic¬ 
itly, let us decompose Hu into an unperturbed part Hq 
and a perturbation V : 

Ho = -dzihdz) + d!if2d!), (E21) 

1 / = -giVi + 52Vl - idzhi)Vl + Vi {/i2,. 

(E22) 

The lowest eigenvalue state for the unperturbed part can 
be expanded as 

uix) = + kzxi + klx 2 + •••], (E23) 

where xi satisfies (E19). Since /i and /2 are invariant 
under z —>■ L — z, the solution of (E19), Xi(z), is an 
odd function, xi(T — z) = —yi(z). Due to the periodic 
boundary condition, Xi(0) = Xi(T) = 0. 


The eigenvalue for the unperturbed part is obtained as 
Eo{kz) = r 1^ (f uHqu = X2*kl + 0{kl). (E24) 

i \w J 


Since V is of order Vi, i.e., k]_ in momentum space, the 
first-order correction SHq gives the contribution of order 
k]_, which is given by the expectation value of V for u: 


5Eo = 


f \u\‘ 


uVu 


kzk]_ ^ |5i(Xi + Xi) - ^2(Xi + xi')} + E){klk\) . 

(E25) 


In the second line, we used (ElO) and integration by 
parts. The integrand is odd under z —>■ L — z be¬ 
cause xi(-^ - z) = -Xi{z), h{E - z) = 5i(z) and 
/i 2 (T — z) = h 2 {z). Therefore, the integral in the sec¬ 
ond line of (E25) vanishes, which proves the absence of 
the kzk^ term in Eo- 
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